Tunable Luttinger liquid physics in biased bilayer graphene 
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Electronically gated bilayer graphene behaves as a tunable gap semiconductor under a uniform 
interlayer bias V g . Imposing a spatially varying bias, which changes polarity from — V g to +V g , 
leads to one dimensional (ID) chiral modes localized along the domain wall of the bias. Due to the 
broad transverse spread of their low-energy wavefunctions, we find that the dominant interaction 
between these ID electrons is the forward scattering part of the Coulomb repulsion. Incorporating 
these interactions and the gate voltage dependence of the dispersion and wavefunctions, we find 
that these ID modes behave as a strongly interacting Tomonaga-Luttinger liquid with three distinct 
mode velocities and a bias dependent Luttinger parameter, and discuss its experimental signatures. 



Graphene has been the focus of intense research in re- 
cent years [l[ due to the rich physics of its massless Dirac 
fermions. Theoretical proposals have shown that spa- 
tially modulating an applied gate voltage in monolayer 
graphene leads to tunable anisotropic transport proper- 
ties For bilayer graphene (BLG), it was experimen- 
tally demonstrated that applying an electric field per- 
pendicular to the layers, by sandwiching it between two 
gates, induces a gap in its electronic spectrum that is 
of the order of the applied interlayer bias 0-0] • Biased 
BLG thus behaves as a semiconductor with a tunable 
band gap. Such tunability of electronic properties opens 
up possibilities for graphene based devices [1 21 ] - 

Recently, Martin, Blanter and Morpurgo Q have pro- 
posed a setup where two neighbouring regions of BLG 
are biased with opposite parity, shown schematically in 
Fig. (jla[) . In this geometry, gapless one dimensional 
(ID) modes were shown to emerge at the interface where 
the bias reverses sign (see also Ref.[H). These modes 
are analogous to domain wall fermions studied in the 
context ofpolyacetylene [§], charge density waves in 
graphene [10j, and field theories in high-energy physics 
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These modes in BLG may also be viewed as switch- 
able nanowires which can be turned on and off using dif- 
ferent gate voltages for a given gate configuration. 

In this Letter, we study the effects of interactions on 
the low-energy modes of a single such wire in the orig- 
inal setup proposed in Ref.pJ. We find that the low- 
energy wavefunctions of the ID modes have a broad 
spread in the direction transverse to the interface — this 
leads to the dominance of the forward scattering part of 
the Coulomb interaction between electrons, in a manner 
akin to large radius carbon nanotubes 13J. Within an 



abelian bosonization framework, incorporating these for- 
ward scattering terms is shown to lead to a strongly inter- 
acting Tomonaga-Luttinger liquid [3]. Remarkably, we 
find the Luttinger parameter of this liquid is tunable by 
adjusting the gate potential. This results from two com- 
peting effects: (i) An increased bias causes further con- 
finement of the wavefunctions to the interface, enhancing 



the effect of interactions; (ii) An increase in the bias in- 
creases the Fermi velocity of the low-energy modes, sup- 
pressing the effect of interactions relative to the kinetic 
energy. The net result is that the Luttinger parameter 
in the total charge channel, K c +, can be varied between 
0.15-0.2 by increasing the bias over an experimentally ac- 
cessible range. At the same time the Luttinger parameter 
in the transverse charge channel, K c - w 0.63, is relatively 
independent of the bias. We thus show that gated BLG 
can realize a tunable Tomonaga-Luttinger liquid. Such 
band structure and wavefunction tuning of Luttinger liq- 
uids has been suggested in a few other systems recently 
— in cold atomic gases 15, 16} . in magnetic waveguides 



in graphene 17(, in carbon nanotubes in crossed elec- 
tric and magnetic fields [HI, [l9[ , and in gated topological 
insulators [20(. Such a Luttinger liquid with dominant 
forward scattering is also expected to arise at charge den- 
sity wave (CDW) domain walls in single layer graphene, 
where the CDW involves a weak sublattice density mod- 
ulation induced by an appropriate substrate fioj . 

Non-Interacting System. - BLG consists of two 
Bernal stacked graphene layers as depicted in Fig. [lb] 
We label the carbon atoms in the bilayer by a unit cell 
index i, a sub- lattice index s = a,b, and a layer index 
£ = 1,2 labelling top and bottom layers respectively. The 
distance between neighboring carbon atoms in the same 




FIG. 1: (a) Schematic diagram of external gates. Vt + Vb 
controls the bilayer doping while Vt — Vb controls the gap 
via the depicted perpendicular electric field, (b) Structure of 
bilayer graphene with a 'bias kink' in V g = Vt — Vb at y — 0. 
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FIG. 2: Dispersion about the K-points with (a) V g = 0.02i 
and (b) V g = 0.08i. Edge-mode bands are indicated by the 
labelled arrows and bulk-states by the hatched region. 




-100 " 100 -100 ~ ~ 100 

y y 

(a) (b) 



FIG. 3: The modulus square of the zero-energy wavefunction 
of (a) the 0-band and (b) the 7r-band on the dominate sites 
a,2 and b\ at V g = 0.02t and V g = 0.08i (independent of being 
a right or left mover). 



layer and on the same sublattice is d ss 2.46A, while the 
interlayer distance is d± » 3. 34 A. The minimal tight- 
binding model for electrons in BLG consists of a nearest- 
neighbor hopping amplitude t « 3eV within each layer, 
and an interlayer hopping amplitude tj_ ss 0.15eV be- 
tween sites (i, s = a,£=l) and (i, s = b,£=2). Henceforth, 
we use units where h=t=d= 1 , and set i = ma + nh 
where a, b are unit vectors depicted in Fig. Ilbl 

In order to describe the effect of the external gates, 
we add a potential term, -|2r i<t (-1) V g (y niSt t)ns. l(T , 
to the tight-binding Hamiltonian, where nR, ff is the elec- 
tron number operator for spin-er at the site labeled by the 
set of indices R = (m, n, s, £). Here, we have assumed 
that the bias V g (y) depends only on the y-coordinatc, 
which is determined by (n,s,£), and is m-independent 
so that translational symmetry is preserved along a 21 1 . 
We assume Vt + Vb (see Fig. 1) is chosen to fix the bi- 
layer chemical potential at the mid-gap. For later con- 
venience, we set R = (to, r), where r = (n,s,£). For a 
general potential profile with V g (y > 0) = — V g (y < 0) 
and V g (y —5- ±oo) = ±V g , the bulk region far from y = 
has a gap A m V g , while the interface has been shown 
to support gapless ID modes Q. Here, for simplicity, we 
assume the gates induce a potential with a step profile 
V g {y) = VgSign(y), since the dominant effect of a poten- 
tial can be shown to come from the str eng th of the bias 
and not the details of its spatial profile [22j . This Hamil- 
tonian can be solved by Fourier transforming along the 
direction parallel to the interface, with momentum la- 
belled k, and then numerically diagonalizing the matrix 
for each k. Our numerical results for various gate biases 
that have been shown to be experimentally viable Q, 
(see Fig. |(5J)) show two right (left) movers for each spin 
at the K (-K) point, consistent with Ref. hi. Drawing an 
analogy with two leg Hubbard ladders [23H25I ]. we refer 
to the high/low energy bands as the 7r/0-bands. 

As seen qualitatively from the figure, and explicitly in 
Fig. ((SJ), the Fermi velocity of the two bands are equal 
and they both change significantly with the bias voltage. 
Fig. |3]) shows the plots of the modulus square of the zero- 
energy wavefunctions in the 0-band and 7r-band at high 
and low gate voltage. We see from here that the spread 



of the wavefunction transverse to the interface direction 
varies significantly. Note that the zero-energy wavefunc- 
tion of the 7r-band is related by interchanging the layers 
and reflecting about the interface so that wavefunction 
leans in to the region of high potential instead. Both, 
this fact and the equality of the Fermi velocities of the 
two bands {Vp = Vf), are a consequence of the sym- 
metry that relates the two band by an inversion about 
the K (-K) point at low energies pf . 

Effective ID Hamiltonian. — To derive the effective 
low-energy ID hamiltonian, we assume a suitable energy 
cutoff that is smaller than the bulk gap and focus on those 
single particle states that lie within this energy window 
and are confined to the ID interface region. To do this, 
we first expand the field operators in the complete ba- 
sis, #r, ct = ^ E fc ,« e ikmd <p%(r) fa, where <(r) is the 
wavefunction of the state in band-a with momentum k. 
We then restrict the bands to the set a = {0,7r} and 
consider only momenta in the vicinity of the four Fermi 
points, ±k F and ±k F . An additional simplification is 
made by neglecting the small momentum dependence of 
the wavefunctions since c/?^ fcQ+9 (r) w (p^ L (r) for small 
momenta q, where <p^^ L (r) is the zero-energy wavefunc- 
tion at ik F . In doing so, the r-dependence of the wave- 
function can be separated to yield the low-energy field 
operators projected to the ID subspace via 



r— ±,a— {0,7r} 



tf(r) e ***"V£r(*). 



(1) 



where r is the label R/L for left/right movers and takes 
the values +/— in the expression, and ip^^x) are slowly 
varying field operators exclusively dependent on the po- 
sition along the interface (which we now denote by the 
continuous variable x = md). 

We now rewrite the entire Hamiltonian in terms of op- 
erators in the reduced ID subspace. The free part is sim- 
ply linearized to give E| g |<A Er„„ rqV F ^(q)fa{q). 
The effective interaction between fermions in the 
ID channel is obtained by a straightforward sub- 
stitution of Eqn. ([1} into the Coulomb term, 



\ E CTCT < £rr< *l,A,*>U(R, R'^rv^r,., 



fol- 
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lowed by a summation over r. This gives rise to various 
scattering terms, many of which are rapidly oscillating 
and can be dropped. The effective Hamiltonian obtained 
contains many terms of the general form 



V 



(i) 



( 2 ) 
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Here, V^ s = v£g yS (rikf - r 4 k 5 F ) is the Fourier compo- 
nent of the effective ID potential 

^? 75 (^-^o=E^ R ' R o<*( r )<V)^ 3 ( r/ )<(r). 

r,r' 

The effective interaction Hamiltonian contains all terms 
of the form in Eqn. © that have a combination of R/L 
and band indices that conserve (crystal) momentum. The 
index i classifies the scattering processes using standard 
g-ology notation [23j-|28| (see Fig. Q): i = 1 refers 
to backscattering, i — 2 to forward-scattering involv- 
ing both right and left movers, and i = 4 to forward- 
scattering involving only right or only left movers. The 
number of distinct processes is greatly reduced by the 
fact that all interband scattering terms with parallel spin 
merely renormalize the coefficients of a corresponding in- 
traband term with parallel spin. 

This Hamiltonian is qualitati vely similar to that ob- 
tained for Hubbard ladders 23|-25| with two significant 
differences. First, both the wavefunctions and Vf are 
sensitive to changes in the applied gate bias. The mod- 
ified wavefunctions alter the distance dependence of the 
effective Coulomb interaction, while the change in Vf ad- 
justs the relative interaction strength parameterized by 
the fine-structure constant a — > ol^-. Second, we note 
that the long-range nature of the Coulomb interactions, 
together with the large spread of the low-energy wave- 
functions, causes the small momentum forward scatter- 
ing processes to dominate. This is reminiscent of large ra- 
dius single wall carbon nanotubes, where the extension of 
the wavefunctions around the tube radius suppresses the 
bare backscattering 13:] . We have checked that the bare 
values of these backscattering and interband scattering 
terms are very small, consistent with this argument. For 
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Vg = 0.02t, so that such processes are expected to be 
important only at very low energy and temperature. We 
therefore first focus on the forward scattering processes. 

Bosonization. — Using the standard abelian bosoniza- 
tion procedure [14| , we introduce the bosonic field 4> aa (x) 
and the phase 9 a<T (x) whose spatial derivative d x aa — 
Tlaa is conjugate to 4> aa (x). Fermion operators can be 
represented in terms of these boson fields via ij)^ a (x) ~ 
e i(r4> atT (x)~§ a „(x)) _ It ig a simp i e ma tter to rewrite the 
density-density interactions in the boson representation 



by means of the relations 9 X <; 
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FIG. 4: Examples of various scattering processes, 
between the bands allows us to lighten our notation by 

defining Va = Vaaaa = Vaaaa and Vb = Vasaa = Vl&aa- 

This leads to the Hamiltonian, 

Ht = ^- fdx(d x <^) T u ■ k-\d x $>) + (d x e) T u ■ K(d x e), 

2ir J 

(3) 

(4) 
(5) 



with 
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Here g A/B = (2V A/B )/V F , and $ 
with a similar definition for 0. 

This Hamiltonian is diagonal in the total/transverse 



basis defined via 



'■>lA) 



± 



where v labels the 



spin (s) or charge (c) sector and <^ c (s)a = <A«t ^ 4>ai- 
In this basis, the spin and charge sectors decouple. The 
spin modes are unaffected by interactions, the Luttingcr 
parameters K s ± = 1 and the velocities u s ± = Vf- The 
charge modes have renormalized velocities and nontrivial 
Luttinger parameters, given by 



u c ± = Vf (1 + y c ±) z 
K c ± = {i + y c ±y~ 2 , 



(6) 
(7) 



d x 9 a 



7!" {pRaa — PLaa)- In addition, the symmetry 



where y c ± — 2(Va± Vb)/ttVf- At the Gaussian level, the 
only effect of the interactions is thus to strongly mod- 
ify K c ±,u c ±. Fig.© shows these parameters plotted for 
various gate voltages. As seen from the figure, K c+ can 
be tuned significantly by the external bias; by contrast, 
A" c _ r; 0.63 (not shown) is relatively bias independent. 

Observable consequences. — The strong interactions 
in the charge channel lead to three different velocities for 
the spin (Vf) and charge (u c ±) modes in the Luttingcr 
liquid. Mapping out these dispersing modes, as has been 
done in semiconductor heterostructures [29j], appears to 
be challenging in the biased BLG system. A more accessi- 
ble signature of the Luttinger liquid physics is the energy 
dependence of the single particle density of states (DOS). 
We expect n(e) ~ e Q , with a > 0. Such a suppression 
of the DOS is expected to lead to a tunneling conduc- 
tance G ~ T a (for voltages eV -C k B T) or a nonlinear 
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scanning tunneling spectroscopy would then be useful to 
uncover the underlying Luttinger liquid physics 31 1. 
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FIG. 5: Fermi velocity Vf and mode velocities u ± in the 
charge sector (in units of td), and the Luttinger parameter of 
total charge sector as a function of V g . We assume a screening 
length of 1000d ~ 0.3fim for the Coulomb interaction and a 
short distance cutoff of 0.5d. 



differential conductance dl/dV ~ V a (for eV > k B T). 
We find a bulk = \{K C+ + + K c _ + - 4) 

and Qfcdgc = + ^c- ~ 2), so that bias depen- 

dent tunneling exponents are expected to be observed. 
Various charge density, spin density and superconduct- 
ing pair correlators are expected to show power-law de- 
cays in this intermediate energy Luttinger liquid regime. 
We find that charge and spin density wave operators 
at 2k F and 2k F are most strongly enhanced by inter- 
actions in this regime, decaying along the interface as 
\x — X '\~( 2+Kc + +Kc -^ 2 , w it n the precise lattice scale mod- 
ulation pattern in (n, s, £) being determined by the pref- 
actors set by the wavefunctions ip^j L (r) from Eq.([T]). 
Modulations at k F ±k F are subdominant. 

Discussion. - - We have shown that BLG in a suit- 
able gate geometry can realize a tunable Luttinger liq- 
uid with four gapless modes. This is expected to break 
down once backscattering and interband scattering terms 
become important; since the bare values of these inter- 
actions are small and that they are all marginal, and 
thus flow slowly, there is an intermediate energy window 
where the physics discussed above should be observable. 
At very low energy (or temperature), however, we ex- 
pect that such processes will gap out all sectors except 
the c+ channel, which should still exhibit Luttinger liq- 
uid physics. Turning to the effect of additional interlayer 
hopping terms (73) between the a2 and b\ sites, we have 
checked that this renders V F ^ V F . This asymmetry 
is small for moderate bias voltages; further, interband 
scattering tends to equalize the velocities 2|| 27 1 so that 
small velocity asymmetries are expected to be unimpor- 
tant. We have assumed that the Fermi level is tuned, via 
Vt + Vb, to be precisely in the middle of the gap - small 
deviations which tend to slightly dope the interface states 
while leaving the bulk gap intact will not qualitatively al- 
ter the physics discussed here. While the 'domain wall' 
modes discussed here are topologically protected Q inde- 
pendent of the precise bias profile transverse to the wire, 
disorder along the wire direction will lead to backscatter- 
ing. Backscattering is somewhat mitigated by the wave- 
function spread, but we expect it will lead to insulat- 
ing behavior at very low energy 30]. Fourier transform 
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